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The energy and baryon densities in heavy ion colhsions are estimated by ana- 
lytically solving a 1+1 dimensional hydrodynamical model with source terms. Par- 
, ticularly, a competition between the energy and baryon sources and the expansion 

O ' of the system is discussed in detail. 

§ ; PACS: 24.10.Nz, 25.75.-q, 12.38.Mh 

^ Theoretical evidence from lattice simulations [jl[ of Quantum chromodynamics (QCD) 

for a phase transition from a hadron gas to a quark-gluon plasma at high temperatures 
has prompted experimental efforts to create this new phase of matter in the laboratory 
^ I during the early stages of relativistic heavy-ion collisions. Besides the deconfinement phase 
^ I transition, the medium effects in hot and dense nuclear matter produced in the collisions 
change the characteristics of strong interaction at hadronic level. For instance, the effects 
^ ■ lead to the transition from the chiral symmetry breaking phase to the phase in which it 
^ . is restored^. Since 1986, oxygen, sulfur, silicon and other beams have been successfully 
accelerated to 14.5 A GeV at the BNL AGS, and to 200 A GeV at the CERN SPS. In 
the past few years heavy-ion collisions are also investigated at the GSI SIS with relative 
low energies of 0.8-1.8 A GeV to study the in-medium effects related to high densities. 
From the above experiments, people have really seen some new phenomena such as the 
J/?/'-suppression|0], the dilepton distribution shifty and the ^ enhancement |5[| . It is, 
however, still difficult to say that definite signals of formation of quark-gluon plasma in 
the collisions are found. This vague situation is expected to be changed in the near future, 
since the BNL RHIC will start to run soon, providing center of mass energies of up to 100 
-|- 100 A GeV and beams as heavy as uranium. To further study the high density effects 
at hadronic level, a new generation of Cooling Storage Ring(CSR) with incident energies 
of 0.4-0.9 A GeV is going to be constructed at LanZhou, China. 

The key point for producing quark-gluon plasma or discussing in-medium effects in 
heavy ion collisions is whether the energy density or baryon density achieved in the col- 
lisions is high enough. From a general consideration, the energy density required for 
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producing the quark-gluon plasma must be several times higher than the energy density 
of normal nuclear matter, and the baryon density for a visible change of chiral properties 
in nuclear matter should be at least the order of the normal nuclear density^. Can such 
a high energy density or baryon density be achieved at SPS and RHIC or at AGS, SIS 
and CSR? The answer to this question depends strongly on the degree of nuclear stopping 
power. The experimental data from AGS|^ demonstrated quantitatively that around and 
below the AGS energy any incident nucleus can be fully stopped in a heavy target. In con- 
trast with this, the data[0 from SPS show that the incident energy is not fully deposited 
in the target, the colliding nuclei are partially transparent. Since the degree of stopping 
power is variable with the incident energy, it is believed that a baryon-rich matter may 
be created either in the central region at lower energies or in the fragmentation region at 
higher energies. 

Assuming full stopping, the maximum energy density and baryon density in the center 
of mass frame for a symmetrical A — A collision can be simply evaluated to be 

emax = 2nNmNCOsh^{yc) , 

rimax = 2nAr COsh(?/c) , (1) 

where n^v is the density of nucleons in normal nuclei, m^v the mass of the nucleon, and yc 
the rapidity of the colliding nuclei in the center of mass frame. 

For the collisions with partial transparency, Bjorken's formula[§] gives 
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where tq is the formation time of the secondary particles, Ra the radius of the colliding 
nuclei, and (^j^ ) the transverse energy distribution in the central region. 

\ "-y / max 

Hydrodynamical models are widely used to describe the macroscopic and collective 
effects in heavy ion collisions, assuming thermal and chemical equilibrium through a large 
number of nucleon- nucleon and secondary collisions. Landau was the first to treat high 
energy collisions with fluid dynamics, and Bjorken^] considered partial transparency in 
relativistic collisions. Assuming the colliding system is boost invariant, the maximum 
energy density reached at time tq decreases monotonously as[§] 

e{r) = tmax (3) 

in the central rapidity region, where Vg is the sound velocity in the fluid. 

In Landau and Bjorken's discussions there is no consideration of the initial collision 
process. Two nuclei are assumed to collide at a unique space-time point (x, t) = (0,0). 
This means that all nucleons in the projectile collide with the target at the same time 
r = 0, and the secondary particles from different collisions are formed at the same time 
r = Tq. Therefore, the maximum energy density emax in (|l]) or (|) is just a simple 
superposition of the N — A collisions. After time tq, the energy density drops down in 
terms of (H) due to the expansion of the system. 
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It is clear that any nucleus-nucleus {A — A) collision with a finite rapidity yc in the 
center of mass frame takes a finite collision time 

d,= ^-^2^^0, (4) 

smh He 

long or short according to the incident energy and the size of the colliding nuclei. Here 
is the nucleon radius. Considering a central ^^'^Au — Au collision at incident energies 
Eiab = 0.4 — 14.5 A GeV, namely from CSR to AGS energies, the collision time 4/m < 
dt < 23 fm is comparable with or even larger than the total time of the fiuid description of 
the system. However, when we consider the same colliding nuclei but at RHIC energy, the 
collision time is only 0.1 fm, the neglect of its infiuence on the fiuid evolution is therefore a 
good approximation in high energy limit. In this paper, we will concentrate our attention 
on the dependence of the energy and baryon densities on the collision time for heavy ion 
collisions with full nuclear stopping. 

The boost invariance in the longitudinal direction used to derive the scaling solution (H) 
is a good assumption at early stages of high-energy nuclear collisions when the transverse 
expansion can be neglected. As the transverse fiow develops with increasing time, the 
fiuid cools down more and more rapidly [^. At the late time in the evolution of the 
reaction, the fiuid will expand radially outward, and one can use boost invariance in the 
radial direction instead of in the longitudinal direction [^. In fact, the radial scaling 
solution is widely used in analytic study of relativistic heavy-ion collisions, for instance, 
the investigation of disoriented chiral condensate Compared with the high-energy 



collisions, the kinematics dominated longitudinal expansion in low-energy collisions is 
much weaker, but the transverse expansion does not change a lot. For example, when 
the colliding energy increases from about 20 A GeV at SPS to 200 A GeV at RHIC, the 



transverse fiow velocity changes from 0.55c to 0.65c|]T2[. On the other hand, what we are 
interested in here is the initial collision period. At high energies, this period is very short, 
the longitudinal scaling solution is a good approximation in the period. At low energies, 
the collision time dt has the same order as the total evolution time of the fiuid. For most 
time of the collision period, the transverse fiow is comparable with the longitudinal fiow. 
Therefore, for analytically estimating the energy density and baryon number density in 
low-energy collisions the radial scaling solution seems better than the longitudinal scaling 
solution. In this case the time evolution of the energy density is replaced by 
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e(T") = emax — ) (5) 



the factor 3 in the exponent arises from the 3-dimensional homogeneous expansion. 

In principle, the expansion of the system starts at the initial time tq. However, in the 
very early stage the expansion is weaker than the homogeneous evolution (^. For the 
sake of simplicity we still use the homogeneous mechanism as an effective description of 
the whole expansion process, but the beginning time of the expansion is shifted from tq 
to Texp- Since the collision time dt is large enough, r^xp can be restricted in the region 

To < Texp < To + dt- 

Let us first analyze qualitatively the behavior of the densities when a non-zero collision 
time is taken into account. We still use N — A collisions as input, but the collision time 
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of each N — A collision is different. This implies a different time in the formation of 
secondary particles. Those produced by the first N — A collision are formed at time tq, 
those by the last one at To + dt, and others between tq and TQ + df. As there is a continuous 
input of secondaries and no expansion before Texp, the energy and baryon densities of the 
system will increase linearly with time in the range tq < r < Texp- When the effective 
expansion begins, there is a competition between the input from the left N — A collisions 
and the expansion of the system. While the former makes the densities increase, the latter 
attenuates the fluid strongly. Therefore, the evolution of the densities in the time interval 
Texp < T < To + (ij depends on the input and expansion rates. After the time tq + dt 
when the input from the N — A collisions ends, the densities drop down monotonously 
like Eq.(|). 

In the absence of viscosity and thermal conductivity, the behavior of a relativistic fluid 
is controlled by the energy-momentum and baryon number conservation lawsl 



df,n'' = a, /i,z/ = 0,1,2,3 (6) 

where T'^'^ = —Pg'^" + (e + P)u^u'^ is the energy-momentum tensor, = nu^ the baryon 
flow, P = P(e) the pressure, e the local energy density and the four- velocity of the 
fluid. The source terms S'^ and a indicate the continuous deposition of energy, momentum 
and baryon number in the fluid due to the N — A collisions. In ultra-relativistic collisions, 
there is almost no net baryon number in the central region, one discusses the energy 
density only. For collisions with full nuclear stopping, we must consider simultaneously 
the energy and baryon densities characterized by (H). 

Using two Lorentz covariant coordinate, the proper time r = y/W—r^ and the space- 
time rapidity ?/ = | In instead of r and t, the above hydrodynamical equations are 
greatly simplifled as 

Ore + 3^—^ = S, 
r 

Ti 

drU + 3— = cr , (7) 
r 

where the source terms on the right side mean the energy density and baryon density 
deposited in the fluid per unit time. Since the total energy input and the total baryon 
number input during the collision time dt are e^axl^^^^ — and nmax{^^^^ — / ^^^^ ■, 
where e^ax and rimax are shown in (P, the source terms are given by 

^ _ e„.ax{A'/' ' I) 
A'lHt 

a = - 



A^l^d 

Taking into account the initial condition 

e(T-o) 
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the solution of the differential equations (0) is 



' e(^o) + (7- - ro)S ro<T <T, 



e r 



exp ) 



^(^0) + - ro)o- To<T < Te^p , 

n(r) = J n{T,,p) (^f + f r (l - (^)') r,,.^ < r < Tq + , (10) 



where the sound velocity Vg comes from the equation of state P = vie. Before the effective 
expansion begins, the energy density increases linearly. In the time period Texp < t < 
To + (ii, the energy density contains two terms. The first one is the time evolution of 
the density at Te^..^, and the second term reflects the contribution of the source term. 
After To + dt there is no more energy input, the second term disappears and the density 
decreases with time. The behavior of the baryon density is similar to that of the energy 
density. Obviously the heavier the colliding nuclei are or the lower the incident energy 
is, the longer is the collision time dt and the more important is the energy and baryon 
diffusion. 

Let us now turn to numerical results. As we have emphasized above, the effect of 
finite collision time on the density evolution is especially important for collisions with full 
stopping. In our numerical calculations below we take ^^"^Au — Au collisions at AGS, 
SIS and GSR energies, namely from 0.4 to 14.5 A GeV. At these low energies, while it 
seems not possible to create quark-gluon plasma in the collisions, the in-medium effects 
at hadronic level are still visible if the energy density or baryon density is larger than 
the relevant density of normal nuclear matter. For instance, from a model independent 
calculation]^ the broken chiral symmetry in the vacuum (T = n = 0) is partially restored 
even in the normal nuclear matter. The energy density and baryon density that are shown 
in the following are scaled, respectively, by the normal energy density = 171^71^ and 
normal baryon density Un = 1/ (^^irr^^, where the nucleon mass itln and the nucleon 
radius are taken to be 0.94 GeV and 1.1 fm. For a perfect fluid, the equation of state 
can be written as P = |e. In real fluid with finite hadron mass, viscosity and thermal 
conductivity, one has P = v^e with t>^ < |. In our calculation we take = |. 

In Fig. la, the energy density, scaled by e^, is shown as a function of the proper time 
r for four values of the expansion time r^xp at the incident energy Eiab = 1.8 A GeV. For 
Texp = To, the effective expansion and the energy input start simultaneously. The energy 
density drops down in the very beginning when the energy input rate is less than the 
expansion rate, and then goes up smoothly when the energy deposited in the fluid is high 
enough. The maximum energy density e(ro + dt)/eN = 0.76 is larger than the initial value 
e(ro)/eAr = 0.67. When the fluid evolves beyond the source region, the energy density 
decreases monotonously. As the effective expansion starts later than the formation of the 
fluid, Texp > To, the energy density increases linearly before the expansion. The increase or 
decrease of the energy density in the region Texp < t < tq + cit depends on the competition 
between the amount of the energy deposited in the fluid and the expansion rate. When 
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Texp is larger than I.Stq, the position of the maximum energy density moves from Tq + dt 
to Texp- For Texp > TQ + dt, there is no expansion in the source region, the maximum energy 
density e(ro + dt) achieves its hmit value (|I|). 

The time evolution of the baryon density, scaled by un, is shown in Fig. lb. The 
structure is very similar to that of the energy density. 

In Fig. 2, we show the maximum energy density and baryon density, scaled by and 
riTv, as functions of incident energy Eiab for the four values of the expansion time r^xp- In 
the case of no expansion in the source region , the maximum energy density and baryon 
density given by (|I]) increase very rapidly and reach 17.5 and 6 at Eiab = 14.5 A GeV. 
When the expansion in the source region is included, the maximum densities are strongly 
suppressed by the expansion. 

Obviously the energy and baryon densities in heavy ion collisions can not reach the 
values estimated in the case of no expansion in the source region. The attenuation of the 
fluid in the initial collision period must be considered. As is well known, the in-medium 
effects in low energy collisions are mainly related to the baryon number density. When 
the in-medium effects become visible, the baryon density should be at least larger than 
the normal nuclear density, Umax/'^N ^ 1- The experimental data of heavy ion collisions 
at SIS with incident energies Eiab = 0.8 — 1.8 A GeV, such as the K~ /K^ enhancement , 
do show the importance of the in-medium effects. This means that the baryon density 
at SIS energies is already beyond the normal nuclear density. From the comparison of 
this analysis with our calculations shown in Figs.l and 2, we can choose approximately 
'To + dt/2 as the beginning time of the effective expansion for the estimation of the energy 
and baryon densities in low energy collisions. 

In summary, we have investigated the influence of non-zero collision time on the evolu- 
tion of energy and baryon densities in heavy ion collisions. When the size of the colliding 
nuclei is large enough and the incident energy is not very high, for instance the colli- 
sions at AGS, SIS and GSR, the initial collision time is long, and the colliding energy 
is deposited in the fluid gradually instead of an instantaneous input in ultra-relativistic 
collisions at RHIG and LHG. In the collision period, there exists a competition between 
the energy and baryon deposition and the expansion. The densities inside and outside the 
source region behavior very differently. They are strongly suppressed by the expansion in 
the source region, and the time to reach the maxima is also remarkably delayed. 

The work was supported in part by the NSFG under grant numbers 19845001 and 
19925519, and by the Major State Basic Research Development Program under contract 
number G2000077407. 
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FIGURE CAPTIONS 

Fig.l 

The time evolution of the energy density (Fig. la) and baryon density (Fig. lb), scaled by 
ejsf and n^v, in ^^''^ Au —^^''^ Au collisions at incident energy Eiab = 1.8 A GeV for four values 
of the expansion time Tg^p. 

Fig.2 

The maximum energy density (Fig. 2a) and baryon density (Fig.2b), scaled by eN and n^, 
as functions of incident energy Eiab for four values of the expansion time Texp- 
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